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Abstract 

Hadron masses are calculated in quenched lattice QCD with a 
fermion action of the D234 type on an anisotropic lattice. Hyper- 
fine splittings for singly charmed baryons are found to be in agree- 
ment with expectations from the quark model and with a magnitude 
slightly larger than experimental values. Masses of doubly charmed 
baryons are also calculated and compared to a variety of model calcu- 
lations. Hyperfine splittings in doubly charmed baryons are found to 
be slightly smaller than in singly charmed states. 



1 Introduction 



Relatively little work has been done on heavy baryons using lattice QCD. 
The most complete study in a relativistic framework is the one done in the 
UKQCD collaboration!]] for both charm and bottom flavoured baryons using 
an 0(a) improved fermion action in quenched approximation. In that work 
hyperfine splittings in both charm and bottom baryons were found to be 
considerably smaller than those predicted by phenomenological models|2|] 
and observed experimentally. A recent NRQCD calculation was able to 
resolve the hyperfine splittings for baryons with b quarks and found values in 
accord with phenomenological expectations. Since the hyperfine splitting is 
an important feature of the baryon spectrum further investigation of charmed 
baryons seems justified. 

In this work we present the results of another quenched lattice simu- 
lation for charmed baryons. Due to limitations in computing resources it 
was not possible to use the same lattice spacing and volume as used in the 
UKQCD calculation. Rather we work on a more coarse lattice (~0.2fm) with 
a highly improved action. Past experience has shown that results of reason- 
able accuracy may be obtained with such lattices 0. In order to check the 
calculation, the spectrum of baryons in the light quark (u,d,s) sector was 
calculated at the same time. As well, meson masses for both heavy and light 
quarks were calculated. The results of all these calculations are in reasonable 
agreement with experimental values and with the results obtained at small 
lattice spacing|| |J. 

The simulation reported here differs from (l| in two other respects. In 
HI the interpolating operators used for the baryons were taken to have a 
form which emphasizes the heavy quark symmetry. Secondly, the spin 1/2 
and spin 3/2 S-like baryons are interpolated by the same operator, a Rarita- 
Schwinger spin-3/2 field. As is well known|0], the correlator of such a field 
propagates both J = 1/2 and J = 3/2 states and it is the J = 1/2 projection 
of this correlator that is identified with the S-like baryons in [[I]]. 

The procedure used here is different. As has been done in the context of a 
QCD sum rule calculation ||, the interpolating fields used for heavy baryons 
are taken to have the same form as those used in the light quark sector. The 
J = 1/2 £ baryon has an interpolating operator which is distinct from that 
used for the J = 3/2 £*, just as the nucleon is usually interpolated by a 
field that is distinct from the A. No assumption is made about heavy quark 
symmetry. An advantage of this is that we can use the same procedure (and 
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computer code) at all masses which provides some check on the results. 

The correlation function used to extract the masses of E's is not obtained 
from the projection of the J = 3/2 field's correlator. In fact, it seems that the 
spin 1/2 projection of this correlator has very poor overlap with the J = 1/2 
ground state. Our finding is that the spin 1/2 projected correlator is very 
small and noisy compared to the correlator calculated directly using a spin 
1/2 field. 

The final conclusion of this study is that in quenched lattice QCD the 
hyperfine splittings for both singly and doubly charmed baryons are in rea- 
sonable agreement with phenomenological expectations. Indications are that 
the splittings may be overestimated compared to experiment which seems 
to be a common tendency for quenched QCD simulations of baryons at all 
quark masses. 

2 Method 

The calculation is done on an anisotropic lattice @ using the gauge field 
action 



S G (U) = (3 



c ps ^(l - -ReTrU ps ) + c rs ^(l - ]-ReTrU r 

ps " rs " 



+c p tJ2( l - ^ReTrU p t) + c rst ^(1 - ]-ReTrU rst ) 

pt ^ rst ^ 



rts & 



(1) 



where ps and rs denote spatial plaquettes and spatial planar 2x1 rectangles 
respectively. The plaquettes lying in the temporal-spatial planes are denoted 
by pt while rectangles with the long side in a spatial(temporal) direction are 
labeled by rst (rts). The c coefficients incorporate the aspect ratio £ = a s /a t 
and gauge link renormalization factors u s and Ut- These renormalization 
factors are estimated using the link expectation value in Landau gauge. 
The fermion action is of the anisotropic D234 type|l(| and has the form 



S F = ^2(cuD u (x) + c 2 iD 2 i{x)) + ^(ci t D lt (a;) + c 2t D 2t (x)) 
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Table 1: Coef ficients appearing in the gauge and fe rmion actions 
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where 

D M (a;) =^(x)(l-^7 i )C^(z)V'(a: + i)+^(a; + «)(l + ^7 i )t/J(a:)^(a;), (3) 

D lt (z) = ^(x)(l - 7 4)f/ 4 (x)^(x + t) + ^{x + t)(l + lA )Ut{x)^{x), (4) 

D 2i (x) = tf(x)(l-fr i )U i (x)U i (x+T)i/>(x + 2t) 

+$(x + 2?)(1 + ^i)U}{x +l)U}(x)^(x) } (5) 

D 2t (x) = ^(x)(l-j A )U4x)U 4 (x + ?)il>(x + 2?) 

+ 2t)(l + 74)^^ + ( 6 ) 

The c coefficients in the fermion action include the aspect ratio, link renor- 
malization and the hopping parameter factors and are shown in Table |l|. 

Hadron masses are calculated from zero-momentum correlation functions 
in the usual way. For mesons the interpolating fields were just the standard 
ones. For baryons some discussion is needed since the procedure used here 
differs from that used in [1]]. Start from the light quark (u,d,s) sector. A 



common choice plfl for the proton operator in terms of u and d quark fields 
is 

e abc [u T a C^d b }u c (7) 
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where a,b,c are colour indices and Dirac indices have been suppressed. For 
A the operator 



-±=e abc {2[«JCV4K + [u T a C^u b ]d c } (8) 

is used. This choice of operators is not unique)]!^] but it allows an easy gen- 
eralization to other baryons [|ll]]. The interpolating operators for the strange 
hyperons S and £* are obtained by the replacement d — ► s in (|7|) and (|8|) 
respectively. Similarly the interpolators for 5 and S* are constructed by the 
replacement u — > s. The A hyperon is more of a problem. In the SU(3) 
flavour limit it would be natural to use the octet lambda 

A 8 = -±=e abc {2[u T a C l5 d b ]s c + [u T a C l5 s b ]d c - [<%C 7B s b ]u c } (9) 

since it is degenerate with the nucleon and A. However, since SU(3) flavour 
is broken, this choice is not compelling. For example, in [|ll| a combination 
A with both SU(3) octet and singlet components was defined. In this work 
a "heavy" A with a form 

[u T a C l5 d h \s c , (10) 

which is natural in the heavy quark limit, is also used. 

The operators used to calculate the masses of the eight ground state 
singly-charmed baryons are taken to have the same structure as the operators 
given above. The E c and S* baryons are obtained by the substitution d — > c 
in (^) and (§) while Q c and Q* are obtained by u — > s, d— > c. For the A c and 
E c the operators (f|) and (|10|) with the replacements s —>■ c and d — > s, s — > c 
are used. Finally, for the remaining two states we take 

K = ^jf hc {[<C l5 c b ]s c + [s T a C l5 c b }u c } (11) 

and 
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^=e abc {[u^C^s b ]c c + [s T a C lllCb }u c + [c T a C^u b ]s c } . (12) 

Rather than using the above operators which have an explicit relativistic 
form one could consider the operators which survive in the limit of a static 
charm quark. The way to do this has been discussed in some detail in the 



context of QCD sum rule calculations|13, 141. However there is no particular 



advantage to taking this limit here. The operators we use contain the leading 



■5 



heavy-quark components so the simulation will decide by itself whether they 
are dominant. The advantage of using the explicit relativistic forms is that it 
allows a unified analysis of hyperfine effects over the whole mass range from 
nucleon and A to charmed baryons. 

The operators such as @ and(|l2|) propagate both spin 1/2 and spin 3/2 
states 0. At zero momentum the correlation function with spatial Lorentz 
indices has the general form 

CiM = (k - llilj)Cs/2(t) + l^C 1/2 (t) (13) 

where the subscripts 3/2 and 1/2 denoted the spin projections. The quantity 
Cs/2(t) was used to extract the mass of the spin 3/2 states. However, it 
was found that the spin 1/2 projection Ci/ 2 (t) was too noisy at large time 
separations to allow for the determination of a mass. 

Hadron correlators were calculated using interpolating operators in local 
form at both source and sink and also applying a gauge invariant smearing to 
the quark propagators at the sink. The Gaussian smearing function, eqn(13) 
of j]i~5|] , was used. Hadron masses were obtained by a simultaneous fit to local 
and sink-smeared correlators. 



3 Results 

The calculations were carried out at /3 = 2.1 on lattices with a bare aspect 
ratio £ of 2. The Landau link tadpole factors were first determined iter- 
atively to be u s = 0.7858 and u t = 0.9472 and these values were used in 
all subsequent calculations. The static potential was determined from both 
spatial and spatial-temporal Wilson loops. From this the lattice spacing and 
renormalized anisotropy were obtained. The results for the lattice spacings 
are a; 1 = (0.977 ± 0.003)GW and aT/ 1 = (1.914 ± 0M7)GeV with a sys- 
tematic uncertainty of O.OlGeV coming from uncertainty in the choice of 
parametrization of the short distance part of the potential [[U|. The renor- 
malized anisotropy was found to be 1.95±0.02 which is compatible with other 



studies done with improved gluon actions at similar lattice spacings|T7| . 

Fermion propagators were calculated on a 10 3 x 30 lattice with Dirichlet 
boundary conditions on the fermion fields. A total of 420 configurations were 
analyzed. With some preliminary tuning it was found that K = 0.182 and 
k = 0.237 gave good values for the J ftp and meson masses so these were the 
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Table 2: Hadron masses for light quarks. Masses are given in GeV, mass 
differences are in MeV. The experimental values in this and other tables are 



k values adopted for the charm and strange quarks for all calculations. Where 
necessary, masses were extrapolated to the physical up and down quark region 
using results from the set of hopping parameter values 0.229, 0.233, 0.237 and 
0.241. The value of the critical k is 0.2429(2). 

First consider the light quark (u,d,s) sector. The pion and p-meson masses 
were fixed at 0.140GeV and 0.770GeV which determines the hopping param- 
eter for up and down quarks (taken to be degenerate) and the lattice scale a p . 
It was found that a~ 1 = (1.99 ± 0.12)GeV which is slightly larger than aj 1 
found from the static potential. This is an inevitable result of the quenched 
approximation. The p-meson mass scale was used in all subsequent calcula- 
tions. The results in the light quark sector are given in Table || with statistical 
errors obtained by a bootstrap procedure. The dominant systematic error, a 
6% uncertainty in the scale determination, is not shown explicitly in this and 
subsequent tables but should be kept in mind. For comparison, results from 
recent calculations (Table II in f| and Table XVII in ||) done at small lat- 
tice spacing and extrapolated to the continuum are also shown. The results 
of our coarse lattice simulation are seen to be quite reasonable. 

In Table 0the A mass calculated with the operator (|9|) is given. The mass 
obtained using (|10|) was essentially identical. This was found to be true for 
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Table 3: Masses for mesons with charm quarks. Masses are given in GeV, 
mass differences are in MeV. 





This work 


Experiment 


Vc 


3.012(4) 


2.980 


J/1l> 


3.087(4) 


3.097 


D 


1.875(6) 


1.867 


D* 


2.007(9) 


2.008 


D s 


1.965(4) 


1.969 


D* s 


2.086(6) 


2.112 


J /if) - Vc 


78(2) 


117 


D* — D 


136(7) 


141 


D* - D s 


122(4) 


143 



all quark masses. 

As check on how well charmed quarks are being simulated we first show 
the results for charmonium and D-mesons in Table || The hyperfine splitting 
between r\ c and J /if) is 78(2)(5)MeK which is in very good agreement with 
the results obtained in jOJ using a completely different fermion action. It 
is considerably smaller than the experimental value which is a well known 
feature of quenched QCD simulations for quarkoniumpO|. The D* — D split- 



tings are compatible with results from NRQCD on similar lattices |21fl and are 
also smaller than experimental values although the suppression of hyperfine 
effects is less pronounced than in charmonium. 

As mentioned in Section 2 the interpolating operators used for charmed 
baryons are taken to have the same form as those used for the light baryons. 
For example, the correlator for S c is calculated directly using operators which 
are the same as used for the strange E hyperon except the mass is increased 
to charm. An alternative is to extract the masses of S-like baryons (e.g., 
S c , VL C ) from the spin-1/2 projection of a correlation function of spin 3/2 
fields( see (|i~3|)). Figure 1 illustrates why this alternative is not used here. 
The correlation functions for fl c and the spin-3/2 projected Q* are plotted 
as a function of lattice time. Also shown is the spin-1/2 projected correlator. 
This spin-1/2 projected piece has a very fast pre-asymptotic falloff and is 
therefore small and noisy in the time region in which one would want to 
extract the mass. This shows that the overlap of the spin-1/2 projection is 
small. It is also worth noting from Fig. 1 that even without any analysis one 
sees that the fl c correlator has a less rapid falloff than that of Q* i.e., Q c is 
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Table 4: Masses of singly charmed baryons. Masses are given in GeV, mass 
differences are in MeV. The experimental values are taken from |TJ| except 
for E' c which is from 





This work 


UKQCD|T| 


Experiment 


A c 


2.304(39) 


2.27 


+3 
-3 


2.285(1) 


K 


2.465(22) 


2.46 ±1 


+5 
-5 


2.453(1) 




2.557(30) 


2.44 t\ 


+4 
-5 


2.518(2) 


'C 


2.454(21) 


2.41 ±1 


+4 
-4 


2.468(2) 


1—1 c 


2.579(14) 


2.51 ±t 


+6 
-6 


2.575(3) 


1— 'c 


2.672(16) 


2.55 ±1 


+6 
-5 


2.645(2) 




2.664(12) 


2.68 ±i 


+5 
-6 


2.704(4) 


n* c 


2.757(14) 


2.66 + 1 


+6 
-7 






91(25) 


-it ig 


+3 
-2 


65(2) 


' — 'c ' — 'c 


94(13) 


-20 lis 


+2 
-3 


70(4) 




94(10) 


-23 t? 4 


+3 
-2 





Table 5: Masses of doubly charmed baryons in GeV 





1 — 'CC 


* 

1 — ^CC 






This work 


3.598(13) 


3.682(20) 


3.697(10) 


3.775(12) 


Potential model |23f 


3.478 


3.61 






Mass formulae [fJjfl 


3.610(7) 


3.735(17) 


3.804(8) 


3.850(25) 


Bag modelp5[ 


3.511 


3.630 


3.664 


3.764 


HQET[|| 


3.610 


3.680 


3.710 


3.760 


Skyrmion model 2*7| 


3.752 


3.793 


3.934 


3.953 



less massive than 

The results for singly charmed baryons are given in Table |j. Overall 
the results are in reasonable close to the experimental values where they are 
known. 

The masses of doubly charmed baryons were also calculated. There are 
no experimental data but a comparison with a selection of model calculations 
is given in Table |5| (A more complete tabulation from various models may 
be found in [24j].) As might expected, without experimental constraints, the 
model calculations vary over a considerable range. The mass splittings from 
our quenched QCD simulation, 84(13)(5)MeV for E* cc -E cc and 78(7)(5)MeV 
for VL* CC — Q cc , are substantial and lie in the middle of the range covered by the 
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models listed in Table ||. They are only slightly smaller than the hyperfine 
splittings found for singly charmed baryons. 

4 Summary 

Hadron masses were calculated with an improved action on an anisotropic 
lattice with a spatial lattice spacing of about 0.2fm. Comparison with sim- 
ulations done at small lattice spacings and extrapolated to the continuum 
indicate that lattice spacing errors are less than 10%. 

The focus of this study is charmed baryons. For both singly and doubly 
charmed baryons spin splittings were found to be in agreement with expecta- 
tions of quark models and other phenomenological approaches. The splittings 
for singly charmed baryons are somewhat larger than experimental values. 
This is in contrast to the small hyperfine effects for singly charmed baryons 
reported by the UKQCD collaboration!]]. 

The calculations reported here were done on a small lattice at a rela- 
tively large lattice spacing. By doing a unified study for light and heavy 
quark masses and comparing to continuum results where possible we have 
some confidence that the correct qualitative pattern of hyperfine effects in 
charmed baryons has been established for quenched lattice QCD. For a pre- 
cise calculation, finite volume and lattice spacing issues have to be addressed. 
It is hoped that this can be done in the near future. 
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Figure caption 

Fig. 1 Correlation function for the fl c field (squares) and the Q* field ( 
spin 3/2 projection (triangles), spin 1/2 projection (circles)) as a function of 
lattice time. 
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